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ABSTRACT 

We present a comprehensive catalog of ultraviolet (HST /STIS and FUSE) absorbers in the low- 
redshift IGM at z < 0.4. The catalog draws from the extensive literature on IGM absorption, and it 
reconciles discrepancies among previous catalogs through a critical evaluation of all reported absorp- 
tion features in light of new HST /COS data. We report on 746 H I absorbers down to a rest-frame 
equivalent width of 12 ml over a maximum redshift path length Az = 5.38. We also confirm 111 
VI absorbers, 29 C IV absorbers, and numerous absorption features due to other metal ions. We 
characterize the distribution of absorber line frequency as a function of column density as a power 
law, dhf/dz oc N-P, where /3 = 2.08 ± 0.12 for VI and /3 = 1.68 ± 0.03 for H I. Utilizing a more 
sophisticated accounting technique than past work, the catalog accounts for ~ 43% of the baryons: 
24 ± 2% in the photoionized Lya forest and 19 ± 2% in the WHIM as traced by VI. We discuss the 
large systematic effects of various assumed metallicities and ionization states on these calculations, 
and we implement recent simulation results in our estimates. 

Subject headings: cosmological parameters — cosmology: observations — intergalactic medium - 
quasars: absorption lines — ultraviolet: general 



1. INTRODUCTION 

The successful connection of recent advances in cosmol- 
ogy to the processes that shape the low-redshift universe 
requires a detailed understandin g of the intergalactic 
medium (IGM). Simulations (e.g..|Cen fc Ostrikerlll999l : 
IDave et al.ll2001t iSmith et al.ll2011[) predict that a large 
portion of the baryonic mass of the universe resides in 
the non-luminous IGM in the form of primordial gas and 
gas processed by galaxies. The systematic characteriza- 
tion of these systems in absorption at different redshifts 
is thus a key step toward understanding the evolution of 
galaxies and the universe in general. 

At high-redshift, the diffuse photoionized IGM, com- 
monly known as the Lya forest, accounts for nearly 
all of the baryonic matter, but simulations and ob- 
servations both find th at this fraction i s reduced to 
~ 30 % by z ~ (e.g., IDave et all 1200 U iPenton et all 
2000). Many of the remaining baryons reside in the 
"warm-hot" IGM (WHIM) phase, which is heated to 
10 5 ~ 7 K as gas is shocked during gravita tional infall, 
cloud -cloud collisions, and by galaxy winds (|Cen fc Fangj 
2006). Because of their low density, these phases of the 
IGM are most readily observed in absorption against 
a background continuum source such as an AGN, us- 
ing Lyman series H I lines as well as highly ionized 
metal species such as C IV and O VI. As these fea- 
tures lie in the rest-frame far-ultravio let (FUV) , re- 
cent l ow redshift (z 5, 0-4) studies (e.g.. IDanforth et all 
2001 iDlmforth fc Shulll [20051 120081 iLehner et all 120071: 



Thorn & Chenl 120081 : iTripp et al J 12008( 1 have focused on 



trograph (STIS) on the Hubble Space Telescope (HST). 
Together, these studies suggest that the column density 
distribution of H I follows a power law, djV /dN cx , 
where j3 ~ 1.7, and that highly ionized metals such as 
O VI follow a similar distribution with (3^2. 

With the installation of the Cosmic Origins Spectro- 
graph (COS) on HST, our sensitivity in the FUV is at an 
all-time high, and we have the opportunity to probe the 
IGM more deeply than ever before. With more than ten 
times the sensitivity of STIS , COS ((Green et all 1201211 
promises to probe the IGM with a larger and more sys- 
tematic selection of sight lines down to lower column den- 
sities. Alre ady, a large number of targets have been ob- 
served (e.g., IDanforth et alll2010 l l2011l ; lNaravanan et all 
120111 iSavage et al.ll2011a|l . A systematic understanding 
of existing STIS and FUSE data in light of the early COS 
data, however, is essential to guiding future IGM studies. 
To this end, we present in this work a systematic cata- 
log of STIS and FUSE IGM absorbers, the largest such 
catalog to date, in light of past studies and new COS 
data. We draw primarily from the extensive literature 
in this field, creating a comprehensive database of IGM 
absorbers that reconciles discrepancies among previous 
lists of absorbers. It contains 746 H I and 111 O VI ab- 
sorbers, as well as detections of many other ion species, 
along 44 AGN sight-lines. We further present new count- 
ing statistics and cosmological parameters derived from 
the catalog, and we compare these results to other work 
in the field. This catalog will lay the groundwork for 
future COS IGM studies. 



observations with the Far Ultraviolet Spectroscopic Ex- 
plorer (FUSE) and the Space Telescope Imaging Spec- 
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2. METHODOLOGY 
2.1. Sight-Line Selection and Data Reduction 

The primary goal of this catalog is to create the 
most comprehensive low-redshift IGM census possible 
with legacy STIS and FUSE data. Therefore, we in- 
cluded all 44 sight lines with zagn ;$ 0.5 for which 
STIS/E140M and/or FUSE data were available at rea- 
sonable quality. The imposed high-redshift cutoff guar- 
antees STIS wavelength coverage of nearly all potential 
Lya features while avoiding complications in interpre- 
tation due to high densities of IGM absorption features 
present in higher-redshift sight lines (z ~ 1 — 2). All data 
have a resolution of 20 km s _1 or better at a signal-to- 
noise ratio per resolution element (S/N) rcs > 5 in most 
cases. The FUSE data provide wavelength coverage over 
905 - 1187 A, while, in most cases, STIS/E140M cov- 
ers 1162 - 1729 A. A few of the sight-lines (3C 273, 
Akn 564, Mrk 509, Mrk 205, and PG 1116+215) lack 
the longest wavelength order of the STIS/E140M data 
and truncate at 1709 A. The data were supplemented by 
COS data in a few cases to assist in the confirmation and 
interpretation of absorption features as described in Sec- 
tion 12.21 These COS data add additional coverage over 
1135 — 1770 A at a resolution of ~ 15 km s _1 and a S/N 
that greatly improves upon the FUSE and STIS data in 
most cases. We used all COS data available for these 
sight lines as of 2011 September 28. Table Q] summarizes 
the final sample of sight lines and data, which cover a 
total Lya redshift pathlength of Az = 5.38. 

The STIS data were reduced as described in our 
group 's previous paper on this subject (jDanforth fc ShulU 
120081 hereafter DS08), and the FUSE d ata were reduced 
as described bv iDanforth et all (|2006l ). We relied pri- 
marily on the hig hest throughput FUSE channel, LiFl 
(jMoos et al J 12000 ). but we also used the LiF2 and SiC 
channels to aid in interpretation and to provide cover- 
age where the LiFl data were absent, suspect, or of 
low S/N. The COS d ata were reduced as described in 
IDanforth etaD (pOlOl ). 

2.2. Catalog Construction 

The creation of the catalog was designed as a critical 
analysis of previously published lists of absorption fea- 
tures in the low-redshift IGM. We began by correlating 
the features reported in a variety of sources, summarized 
in Table [2j for all sight lines in our sample. Although 
some of the sight lines were analyzed in only one previ- 
ous work, many have been studied by two or more groups. 
One sight line, Ton S210, has data from FUSE, STIS, and 
COS, but it has not been included in any major, multi- 
ion IGM absorption survey. As such, our analysis of that 
sight line is original. Each feature was individually in- 
spected in the data, and its interpretation, in terms of 
both structure and identification, was confirmed in con- 
text with related lines. We used any available COS data 
to aid in the interpretation of each feature, but the mea- 
surements themselves were not made from the COS data. 
Although the COS data sometimes refute the existence 
of features that are seen by the other instruments, the 
features were nevertheless included in the catalog so as 
to maintain a homogeneous sample; they were, however, 
flagged as suspect. For cases in which our interpretation 



and all existing literature measurements agreed in the ba- 
sic interpretation of the feature, we computed consensus 
values for the observed wavelength, rest-frame equivalent 
width and Doppler 6-values as means of the measure- 
ments reported in the literature. One-sigma errors were 
computed assuming that the measurements were inde- 
pendent. This assumption, while not formally accurate, 
is roughly approximated by differences in data reduction, 
continuum placement, and measurement technique. As a 
complementary measure of the error, we also computed 
the mean absolute deviation among the individual mea- 
surements. In cases where one or more of the sources 
disagreed on the interpretation of the feature, we criti- 
cally evaluated the arguments and data and either chose 
the best interpretation or, if necessary, remeasured the 
feature before proceeding as in the case of total agree- 
ment. If appropriate, alternate identifications and struc- 
tural interpretations were noted. Throughout this pro- 
cess, a number of previously unreported features were 
found scrcndipitously. Although we measured them and 
included them in the catalog, we note that this project 
was not a systematic search for new absorption systems, 
owing to its literature-based nature. Every redshift with 
an H I absorption system, however, was checked for cor- 
responding absorption in each of the metal ions listed in 
Tabled 

When we measured a feature's properties, we used 
either the apparent op tical depth (AOD) method 
(jSembach fe Savage 1992) or a Voigt profile fit (see dis- 
cussion in IDanforth et all 120061 and sources therein). 
We qualitatively judged which method was most ap- 
propriate for each feature, taking into account the two 
methods' relative strengths. Although the two methods 
usually generate comparable results, the AOD method 
performs better for weak lines or in extremely noisy 
regions of the spectra, while profile fitting is prefer- 
able for strong or blended l ines o r lines with compli- 
cated component structure (|DS08l ). Where n ecessary, 
we ad opted the atomic parameters given by iMortonl 
(2003 ) or, for lines not included in that work, from th e 
INIST Atomic Spectra Database! (fRalchenko et al.112010) . 

We took a conservative approach to the splitting of 
absorption features into separate components. In cases 
where clear separation was visible or multiple lines of 
the same absorber showed consistent asymmetry, we fa- 
vored multiple-component interpretations. On the other 
hand, if separation into components was ambiguous, we 
favored a single-component interpretation, even if the 
overall profile shape suggested the presence of multiple 
components. Where possible, we used high-S/N COS 
data to further inform our interpretation of the compo- 
nent structure. This approach leaves the catalog with 
few, if any, false positives in terms of measured com- 
ponents, but some features may have more components 
than reported here. For comparison, the features re- 
ported in this catalog yi eld cou nting statistics inter- 
mediate between those of DS08L which generally used 
single-co mponent interpreta tions, and the profile-fitting 
results of iTripp et al.l ([2008), who generally used many- 
component interpreta tions. Our componen t splitting is 
comparable to that of iThom fc C hen (2008). Differences 
in component-splitting approaches are expected to have 
little impact on the calculation of baryon densities be- 
cause the primary methods used to derive those quanti- 
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TABLE 1 

44 AGN SlGHTLINES AND DATA USED IN THE CATALOG 



Sight Line 




R.A. 


Decl. 


^AGN 


STIS 


FUSE 


COS/G130M 


COS/G160M 




(J2000.0) 


(J2000.0) 




(ksec) 


(ksec) 


(ksec) 


(ksec) 


Mrk335 


00 


06 


19.5 


+20 12 10 


0.0258 


17.1 


97.0 


2.4 


1.6 


IZwl 


00 


53 


34.9 


+12 41 36.0 


0.0607 




38.6 






TonS180 


00 


57 


20.0 


-22 22 59.3 


0.0620 




16.6 






TonS210 


01 


21 


51.5 


-28 20 57 


0.1160 


22.5 


53.3 


5.0 


5.5 


Fairall 9 


01 


23 


46.0 


-58 48 23.8 


0.0461 




38.9 






HE 0226-4410 


02 


28 


15.2 


-40 57 16 


0.4950 


43.8 


33.2 


6.7 


7.8 


NGC985 


02 


34 


37.8 


-08 47 15.6 


0.0431 




68.0 






PKS 0312-770 


03 


11 


55.2 


-76 51 51 


0.2230 


8.4 


5.5 






PKS 0405-123 


04 


07 


48.4 


-12 11 37 


0.5726 


27.2 


71.1 


22.2 


11.1 


Akn 120 


05 


16 


11.4 


-00 08 59.4 


0.0331 




56.2 






HS 0624+6907 


06 


30 


02.5 


+69 05 04 


0.3700 


62.0 


112.3 






VII ZwllS 


07 


07 


13.1 


+64 35 58.8 


0.0797 




198.6 






PG 0804+761 


08 


10 


58.5 


+76 02 41.9 


0.1000 




174.0 


6.7 


6.3 


Ton 951 


08 


47 


42.5 


+34 45 03.5 


0.0640 




31.9 






PG 0953+414 


09 


56 


52.4 


+41 15 22 


0.2341 


8.0 


72.1 






Ton 28 


10 


01 


02.5 


+28 55 35 


0.3297 


33.0 


11.2 






3C 249.1 


11 


04 


13.7 


+76 58 58 


0.3115 


68.8 


216.8 






Mrk421 


11 


01 


27.3 


+38 12 32.0 


0.0300 




83.9 


1.7 


2.1 


PG 1116+215 


11 


19 


08.6 


+21 19 18 


0.1765 


26.5 


77.0 






PG 1211+143 


12 


14 


17.7 


+14 03 13 


0.0809 


42.5 


52.3 






PG 1216+069 


12 


19 


20.9 


+06 38 38 


0.3313 


5.8 


12.0 






Mrk205 


12 


21 


44.0 


+75 18 38 


0.0708 


62.1 


203.6 






3C 273 


12 


29 


06.7 


+02 03 09 


0.1583 


18.7 


42.3 






Q 1230+0115 


12 


30 


50.0 


+01 15 23 


0.1170 


9.8 


4.0 


11.1 


8.1 


PG 1259+593 


13 


01 


12.9 


+59 02 07 


0.4778 


95.8 


668.3 


9.2 


11.2 


PKS 1302-102 


13 


05 


33.0 


-10 33 19 


0.2784 


4.8 


142.7 


7.0 


6.9 


Mrk279 


13 


53 


03.4 


+69 18 30 


0.0305 


54.6 


228.5 


2.2 


2.7 


NGC 5548 


14 


17 


59.5 


+25 08 12 


0.0172 


69.8 


55.0 


1.9 


2. i 


Mrk 1383 


14 


29 


06.6 


+01 17 06 


0.0865 


10.5 


63.5 






Mrk817 


14 


36 


22.1 


+58 47 39.5 


0.0313 




189.9 


3.4 


3.0 


Mrk 478 


14 


42 


07.5 


+35 26 22.9 


0.0791 




14.2 






PG 1444+407 


14 


46 


45.9 


+40 35 06 


0.2673 


48.6 


10.0 






Mrk 290 


15 


35 


52.4 


+57 54 09.3 


0.0296 




12.8 


3.9 


4.8 


Mrk 876 


16 


13 


57.2 


+65 43 10 


0.1290 


29.2 


46.0 


12.6 


11.8 


3C 351 


17 


01 


41.4 


+60 44 31 


0.3719 


77.0 


141.9 






H 1821+643 


18 


21 


57.3 


+64 20 36 


0.2970 


50.9 


132.3 


0.6 


0.5 


PKS 2005-489 


20 


09 


25.4 


-48 49 53.9 


0.0710 




49.2 


2.5 


1.9 


Mrk 509 


20 


41 


09.7 


-10 43 25 


0.0344 


7.6 


62.3 


9.0 


16.5 


II Zwl36 


21 


32 


27.8 


+10 08 19.4 


0.0630 




22.7 






PHL1811 


21 


55 


01.5 


-09 22 25 


0.1900 


33.9 


75.0 


3.5 


3.1 


PKS 2155-304 


21 


58 


52.0 


-30 13 32 


0.1160 


10.8 


123.2 






Akn 564 


22 


42 


39.3 


+29 43 31 


0.0247 


10.3 


60.9 


1.7 


2.4 


MR 2251-178 


22 


54 


05.8 


-17 34 55.0 


0.0644 




54.1 






NGC 7469 


23 


03 


15.6 


+08 52 26 


0.0163 


22.8 


44.3 


1.9 


1.8 



TABLE 2 

Sources of data used in the catalog 



Number a 



Reference 



Number a 



Reference 



1 
2 
3 
4 
5 
6 
7 
8 
9 
10 
11 



Willigcr ct al. (2006) 
Danforth fc ShuII (20TO 
Thorn fc C hen (2008) 
Prochaska et al. (2004) 
Lehner et al. (20071 
Trip p et al. (2008 . Table 2) 



ScmbachetaL 
^ripj3^t_aL 
Oceerle_et_al. 

AraciletaL 



(2001) 
2001) 
2000) 
2006) 
2002) 



12 Scmbach ct al. (2004) 

13 Richter et al. (2004) 

14 Jenkins et al. (2003) 

15 Narayanan et al. (2011) 

16 Tripp et al. (2008. Table 3) 

17 Savage et al. (20051 

18 Danforth et al. (2006) 

19 Howk et al. (20091 

20 Richter et al. 

21 Tumlinson ct al. 

22 Lehner et al. 



a Identifying number used in the catalog. 
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TABLE 3 
Primary IGM Diagnostic Lines 



Ion 


A rcst 


J 




£ C 

J ion 




log 1 max 




f A) 

( A ) 




CIE 






HI 


1215.67 


0.4164 


0.00 






< 4.19 d 




1025.72 


0.07914 


0.00 






< 4.19 d 


Fell 


1144.94 


0.0830 


-4.50 


0.79 


4. 


.05 - 4.27 d 




1063.18 


0.0547 


-4.50 


0.79 


4. 


.05 - 4.27 d 


CII 


1334.53 


0.128 


-3.57 


0.97 


4. 


.15 - 4.66 d 




1036 34 


0.118 


-3.57 


0.97 


4, 


.15 — 4.66 d 


Ol 11 




118 
1 . lo 


A AQ 


u.y / 








1193.29 


0.582 


-4.49 


0.97 




< 4.30 d 


SII 


1259.52 


0.0166 


-4.88 


1.00 




< 4.47 d 




1253.81 


0.0109 


-4.88 


1.00 




< 4.47 d 


Felll 


1122.52 


0.0544 


-4.50 


0.89 


4. 


.27-4.58 d 


Si III 


1206.50 


1.63 


-4.49 


0.90 


4. 


.30 - 4.85 d 


cm 


977.02 


0.757 


-3.57 


0.83 


4, 


.66 - 5.01 d 


sin 


1012.50 


0.0438 


-4.88 


0.84 


4 


.47 - 4.93 d 




1190.20 


0.0237 


-4.88 


0.84 


4. 


.47 - 4.93 d 


Si IV 


1393.76 


0.513 


-4.49 


0.35 




4.8 e 




1402.70 


0.254 


-4.49 


0.35 




4.8 e 


CIV 


1548.20 


0.190 


-3.57 


0.29 




5.0 e 




1550.77 


0.0948 


-3.57 


0.29 




5.0 C 


SIV 


1062.66 


0.0494 


-4.88 


0.61 


4. 


.93 - 5.13 d 


NV 


1238.82 


0.156 


-4.17 


0.24 




5.25 e 




1242.80 


0.0777 


-4.17 


0.24 




5.25 e 


OVI 


1031.93 


0.1325 


-3.31 


0.22 




5.45 e 




1037.62 


0.0658 


-3.31 


0.22 




5.45 c 



Absorption lii 



Solar 
Peak 



photosphc 
on fr actio 



cillator strength 
; clement abund; 



2003). 



co, log(X/Jf) fTspT 
CIE lijutherlamOny 



nd ct al. 2009). 



CIE; tej 



CIE lijutnerland fc DopitaHT^iD . 

turc is peak CIE abundance. 



ties depend only on the total column density observed. 

After compiling complete lists of the properties of in- 
dividual absorption lines from a given system, we took 
a similar approach in determining the consensus proper- 
ties of the absorber inferred from those lines. In cases 
where only a single line was present for the absorber, the 
Doppler 6-value was taken to be the same as the mean 
value for the feature. The column density was computed 
as the geometric mean of the literature values, with er- 
rors treated as above. In many cases, more than one 
absorption line was present for the absorber. For metal 
species with multiple lines, we proceeded as with single- 
line detections but used a mean weighted by the expected 
relative strengths of the observed absorption lines. If any 
of the lines was deemed suspect, for example because of 
excessive noise or blending, it was not included in the 
mean. The H I absorbers were treated in a more com- 
plicated fashion, because, in addition to allowing mea- 
surements of the individual Lyman series lines, the ion is 
particularly conduci ve to curve-of-growth (COG) analy- 
ses (see discussion in Danfor th et alJl2006t ) when multiple 
lines are present. In cases of agreement among the var- 
ious literature sources, we were able to compute simple 
mean values for the column density and b- value of the ab- 
sorber as above. However, if any remeasurement or rein- 
terpretation was necessary at the level of individual ab- 
sorption lines, we used either new COG analyses or new 
weighted means of the line measurements to determine 
the Doppler 6-values and column densities for the ab- 
sorber. In these consensus mean at the absorber level, we 
also included measurements from literature sources that 
did not report properties for the individual absorption 



lines. In a few cases with extremely high column density 
H I, we estimated the absorber parameters by fitting the 
integrated absorptin from higher Lyman lines (i.e., the 
Lyman decrement), if sufficiently short- wavelength data 
were available. The vast majority of H I absorbers, how- 
ever, show only Lya absorption, and in those cases the 
absorber-level properties were inferred directly from that 
single line. 

Because the various literature sources differ slightly in 
their approach to reporting their measurements, we made 
minor adjustments t o our general approach in individual 
cases. iTripp etlrtl ()2008l ) report both total equivalent 
width and column density for each absorption complex as 
well as profile-fitting measurements of individual compo- 
nents; we used the relevant measurements as appropriate 
to our interpretation of the s tructure of th e featu re. In 
using the measurements of the lLehner et al.l (|2006[ ) study 
of the HE 0226-4110 sight line, we adopted the AOD 
equivalent width and Doppler b- value measurements and 
the profile-fitt i ng co lumn density measurements. The 
IWilliger et aD (l200l study of the PKS 0405-123 sight 
line employed a methodological approach that resulted 
in a large number of very weak H I detections which, 
from an examination of the COS data, seem to be noise 
features rather than real absorption features. While we 
have included these measurements in the catalog for the 
sake of completeness, their equivalent widths and column 
densities have been flagged as upp er limits to avoid bi - 
asing the sample set. Similarly, theHtichtcr ct ldl ()2004[ ) 
study of the PC 1259+593 sight line reports a large num- 
ber of "probable" Lya detections. All were re-evaluated 
in light of the much higher S/ N COS data and inclu ded 
in our catalog, but, as with the IWilliger et al.l (|2006|) pa- 
per, we flagged those that seem be noise as upper limits 
rather than measurements. Finally, iRichter et al.l (2004) 
report a variety of measurements for each detected fea- 
ture. For metal features, we favored the profile-fitting 
measurements where available and otherwise used the 
AOD measurements. For H I, we used the profile-fitting 
measurements for individual features and, when possible, 
the COG measurements for ion-level absorber properties. 

A number of features in the catalog are reported as 
upper limits, owing to excessive blending, weakness rela- 
tive to noise, or for the reasons described in the previous 
paragraph. In these cases, we reported the most con- 
servative limit from the references, which may be either 
a measured value in the case of blended lines or a limit 
based on the noise characteristics of the data. Because 
upper limits arc treated somewhat differently among dif- 
ferent authors, we report a minimum S/N-based equiva- 
lent width for each feature for the sake of consistency: 



W n 



4A 



7200(S/N). 



-(1 + z)- 1 . 



fl 



This quantity roughly corresponds to a 4 a significance 
level in equivalent width for the data in the rest frame 
of the absorber. It assumes a resolved feature with b — 
25 km s , roughly the median b- value for H I and O VI, 
with the factor of 7200 corresponding to the resolution, 
R = A/AA, required to resolve a line with that 6-value. 
The signal-to-noise, (S/N), was calculated as the ratio 
of mean flux to its standard deviation in the smoothed 
continuum. Features with equivalent widths near or be- 



The Low-z IGM 



5 



low the reported value of Wmi n should be viewed with 
skepticism and treated as upper limits. Note, however, 
that these 4 a estimates do not account for fixed-pattern 
noise, which is mitigated by the smoothing applied dur- 
ing the data reduction, or for other systematic effects. 

Table |4] presents the complete catalog of individual ab- 
sorption lines. Columns (1) through (4) identify the fea- 
ture according to sight line, observed wavelength, red- 
shift, and transition, respectively. Column (5) lists an 
"ion code" that correlates absorption lines within a sight 
line in Table 0] with the associated absorber properties 
in Table \5\ Columns (6) through (9) describe the mea- 
sured strength of the feature. Column (6) is a flag that 
identifies the rest-frame equivalent width listed in Col- 
umn (7) as a normal measurement ( "n" ) , an upper limit 
("<"), a lower limit (">"), or an order-of-magnitude es- 
timate ( " ~" ) when a reliable measurement could not be 
obtained. Column (8) is the mean absolute deviation of 
the measurements used to generate that consensus value. 
Column (9) identifies the sources or techniques used to 
construct the consensus value, with every two-digit por- 
tion of the value in the table corresponding to a source 
listed in Table [2] Other numbers are used to identify 
our own measurements: 23 corresponds to a Voigt fit, 
24 corresponds to an AOD measurement, 25 corresponds 
to a COG measurement, and 26 corresponds to a Ly- 
man decrement estimate. Columns (10) through (13) re- 
port on the velocity width of the feature in terms of the 
Doppler ^-parameter in the same format in which the 
equivalent width measurements were reported. Column 
(14) lists the quantity W m i n as defined in Equation Q] and 
evaluated at the observed wavelength. Column (15) lists 
possible alternate line identifications in the rare case that 
identification was ambiguous. Finally, Column (16) flags 
inconsistencies between the COS data (where available) 
and the STIS and/or FUSE data. No flag is present if 
COS data are absent or in agreement with the other data, 
while flags are listed for features that are entirely absent 
from the COS data ("a"), weaker or partially absent in 
the COS data ("p"), or that display a different shape or 
structure in the COS data ("d"). 

Table [5] presents the complete catalog of absorbers. 
Columns (1) through (3) identify the absorber via the 
sight line, redshift, and ion species, respectively. Col- 
umn (4) reports the ion code that correlates Tables [4] 
and El as described in the the previous paragraph. Us- 
ing a similar format to the measurements in Table 01 
Columns (5) through (7) report on the column density 
of the absorber, with Column (5) identifying the type of 
measurement with the same set of flags used in Table [SI 
Column (6) reporting the column density, and Column 
(7) reporting the references or techniques used to gener- 
ate the measurement, again as in Table 01 Columns (8) 
through (10) report on the measurements of the Doppler 
6-parameter in the same manner. 

2.3. Absorber Statistics 

In com piling absorber statistics, we followed the ap- 
proach of DS08. We focus on fourteen ions, listed in Ta- 
ble[3]along with some of their properties. For all statistics 
and derived quantities, we included only absorbers in the 
redshift range 500 km s _1 < cz < (czagn — 1500 km s -1 ) 
and further required that z < 0.4. These limits help to 
exclude absorbers intrinsic to either the Milky Way or the 



background AGN while also ensuring that the sample is 
restricted to the local, low-redshift universe where sensi- 
tive H I measurements via the Lya transition can be ob- 
tained. We further excluded all absorbers at z > 0.104 in 
the Q1230+115 sight line (zagn = 0.1170 ) because these 
featu res are likely intrinsic to the AGN (Ganguly et al. 
2003j). Despite these velocity limits, some high-velocity 
AGN outflow systems may contribute to our sample. 
These limit s are l ess restrictive than those adopted by 
i Tripp et al.l (|2008f l and comparable to those adopted by 

For each ion species, we constructed simple histograms 
of the number of absorbers in 0.2 dex bins in log A, as 
well as distributed histograms in wh ich each detection is 
characterized by a Gaussian (|DS 08j). We proceeded to 
calculate the number of absorbers per unit redshift, 



cW(logA) = ^M(logA i ; 



dz 



Az, 



(2) 



This quantity is a function of column density, A, be- 
cause it must be corrected for incompleteness. The cor- 
rection is acheived by dividing the number, Mi, per bin 
by an effective redshift path length, Azi, which is itself 
a function of column density and accounts for variations 
in sensitivity to the absorption feature in question. 

We establish this redshift pathlength by considering 
the S/N in the continuum of the highest-S/N data cover- 
ing a given wavelength. This (S/N)* vector was rebinncd 
to a resolution comparable to the data. The smoothed, 
normalized flux was then used to mask (S/N) * such that 
areas of saturated absorption were set to (S/N)* = 
while regions of unsaturated absorption were scaled with 
the optical depth, r*, as (S/N)* oc (r* + l) -1 . This 
approach effectively removes lines of any type from the 
effective pathlength. The resulting vector was converted 
to a 4 a minumum equivalent width vector defined, as in 
DS08, by 



4A 



i?(S/N); 



(3) 



where the instrumental resolution, R = A/ A A, was set 
to 15,000 or 42,000 for FUSE or STIS data, respectively. 
The effective redshift pathlength for singlet lines is then 
simply the total path length over all sightlines such that 
W\ > Wmin- This is easily extended to multiplet lines 
of different line strengths (/A) by requiring, in the rest 
frame of the weaker line, that either W\ > W m i n (zi) or 
W x > [W^ in {z 2 )]{fx\xlhM). 

Uncertainty in dM /dz is dominated by the uncertainty 
in a bsorber numb er (Mi) per bin due to Poisson statis- 
tics (|Gehre ls 1986). Inherent uncertainty in effective red- 
shift pathlength Azi is negligible, but this quantity con- 
tributes uncertainty arising from its variation over the 
width of a column density bin. We approximate this un- 
certainty as 



d(Az l ) = - [{Az) l+1 - (Az)i-i] 



(4) 



as suggested bv lDSOSi in which these two sources of error 
were added in quadrature. 
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TABLE 4 

IGM ABSORPTION LINE PROPERTIES 



Sight 


^obs 


z 


Line 


Ion 


W 


W 


Wmad 


Ref. 


b 


b 


?>MAD 


Ref. 




Alt. 


COS 


Line 


(A) 






# 


Flag 


(mA) 


(mA) 




Flag 


(km/s) 


(km/s) 






ID 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 



3C249 


1596.95 





31364 


Lya 


1 


n 


148 ±8 




06 


n 


13 ± 1 


23 


21.0 


3C249 


1591.56 





30920 


Lya 


2 


n 


30 ±3 




02 


n 


8±1 


02 


29.0 


3C249 


1590.21 





30809 


Lya 


3 


n 


262 ± 17 




02 


u 


33 ± 1 


02 


29.0 


3C249 


1589.95 





30788 


Lya 


4 


n 


68 ±19 




02 


n 


22 ±4 


02 


29.0 


3C249 


1573.39 





29426 


Lya 


5 


n 


40 ±8 




02 


n 


10 ± 1 


02 


28.9 


3C249 


1544.16 





27021 


Lya 


6 


n 


73 ± 16 




02 


~ 


43 


02 


47.9 


3C249 


1539.82 





26664 


Lya 


7 


n 


96 ±8 




02 


n 


67 ± 19 


02 


51.2 


3C249 


1537.49 





26473 


Lya 


8 


n 


210 ± 19 




02 


n 


33 ±2 


02 


44.5 


3C249 


1532.51 





27021 


SUII1207 


9 


u 


18 ±6 




02 


~ 


19 


02 


44.2 


3C249 


1515.64 





24676 


Lya 


10 


n 


452 ± 16 


17 


0206 


n 


39 ± 1 


02 


51.1 


3C249 


1512.24 





24396 


Lya 


11 


n 


473 ± 23 




02 


u 


56 ±2 


02 


55.4 


3C249 


1511.47 





24332 


Lya 


12 


n 


317 ±27 




02 


u 


53 ±3 


02 


55.4 


3C249 


1510.34 





24239 


Lya 


13 


n 


176 ± 20 




02 


n 


41 ±3 


02 


55.4 


3C249 


1505.34 





23828 


Lya 


14 


u 


164 ± 25 




02 


u 


14 ± 1 


02 


56.1 


3C249 


1504.27 





24681 


SUII1207 


15 


u 


22 ±9 




23 


u 


7±3 


23 


55.7 


3C249 


1459.93 





20093 


Lya 


16 


n 


78 ±18 




02 


n 


12 ±2 


02 


57.2 


3C249 


1459.63 





20068 


Lya 


17 


u 


312 ± 11 




02 


u 


34 ±3 


02 


57.2 


3C249 


1458.69 





17371 


NV1243 


18 


< 


17 




02 








58.4 


3C249 


1455.10 





19695 


Lya 


19 


n 


57 ±4 




02 


~ 


43 


02 


42.5 


3C249 


1454.02 





17371 


NV1239 


18 


u 


29 ±4 




02 


u 


24 ± 7 


02 


43.3 


3C249 


1441.77 





18599 


Lya 


20 


n 


53 ±10 




02 


n 


23 ±4 


02 


35.9 


3C249 


1429.54 





17593 


Lya 


21 


n 


238 ± 36 




02 


n 


39 ±6 


02 


49.8 


3C249 


1429.06 





17553 


Lya 


22 


n 


244 ± 39 




02 


n 


41 ±6 


02 


49.8 


3C249 


1426.84 





17371 


Lya 


23 


n 


49 ± 12 




02 


u 


29 ±5 


02 


27.0 


3C249 


1423.54 





17099 


Lya 


24 


n 


34 ± 10 




23 


n 


23 ±5 


23 


27.0 


3C249 


1412.80 





17099 


SiIII1207 


25 


u 


16 ±3 




02 




10 


02 


31.5 


3C249 


1410.87 





16057 


Lya 


26 


n 


73 ± 10 




02 


n 


29 ±3 


02 


31.2 


3C249 


1396.38 





24396 


Fe III1 123 


27 


u 


16 ±3 




02 




4 


02 


24.8 


3C249 


1391.01 





11925 


NV1243 


28 


< 


10 




02 








32.0 


3C249 


1388.09 





14183 


Lya 


29 


n 


278 ± 10 




02 


n 


38 ± 1 


02 


31.3 


3C249 


1386.55 





11925 


NV1239 


28 


u 


11 ± 2 




02 




15 


02 


32.5 


3C249 


1384.73 





13907 


Lya 


30 


u 


27 ±9 




02 


u 


25 ±7 


02 


31.9 


3C249 


1379.42 





13470 


Lya 


31 


n 


64 ±3 




02 


n 


65 ±2 


02 


25.1 


3C249 


1368.50 





12572 


Lya 


32 


u 


373 ± 21 




02 


u 


31 ± 1 


02 


28.7 


3C249 


1368.07 





12536 


Lya 


33 


u 


151 ± 21 




02 


u 


43 ±5 


02 


28.7 


3C249 


1363.05 





31364 


OVI1038 


34 


n 


22 ± 11 




06 


n 


14 ±2 


23 


24.7 



Note. — Tabic 4 is published i 
. complete description of the co 



edition of the Astrophysit 
nd in Section [2*21 of the tc: 



regarding its for: 



2.4. Baryon Fraction 

Studies of the cosmic microwave background 
(jKomatsu et alj 1201 If) have provided accurate mea- 
surements of the cosmological baryon density, ttf, — 
0.0455 ±0.0028, as a fraction of the critical density of the 
universe, p cx = (3H$/8ttG) = (9.205 x lO" 30 g cm- 3 )/if 
for a Hubble constant, Hq = (70 km s _1 Mpc _1 )/i7o- 
The majority of this matter, both at high and low 
rcdshift, must reside in the IGM. A complete accounting 
of IGM baryons at low redshift is therefore best studied 
with a combinati on of FUV and X-ray spectroscopy 
(Shull et al.l 1201 il ). We address the former here with 
our catalog, which traces IGM gas at T < 10 6 K, by 
considering H I absorbers in the Lya forest and metal 
species as tracers of the WHIM. 



2.4.1. Metal Ions 

We present two quantities that measure baryon den- 
sity contributions as traced by metal ions. In both cases, 
we use five equally spaced bins in z, and we assume 
a mean atomic mass , p = 1.32raH, for Yh c = 0.2477 
(jPeimbert et al.|[2007f) . The contribution to the critical 



density by the IGM mass of an ion is given by 



CPcr 



d/V(log N) 
dz 



Nd(\ogN)dz 

la JN mia 
= (1.365 x 10~ 23 cm 2 )/^ 1 (m ion /amu) 

log JV„ 

f i : \ i 1 1 id .■■ v : 



cW(log N) 




(5) 



where mi on is the ion mass and the integration is over 
column density (cm~ 2 ) and redshift path length. We 
report the uncertainty in f2i on as the asymmetric Poisson 
uncertainty in dAf evaluated at the mean dz, (z), and (N) 
weighted by dAf. 

While fiion is the observed contribution to the critical 
density by a given ion, we are also interested in estimat- 
ing the total baryonic density of all species as traced by 

a particular ion. We denote this estimate as it 
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TABLE 5 

IGM ABSORBER PROPERTIES 



Sight Line 


z 


Ion 


Ion code 


N Flag 


\ogN 
(N in cm -2 ) 


Ref. 


b Flag 


b 

(km s — 1 ) 


Ref. 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


3C249 


0.31364 


H 1 


1 


n 


13.976 ± 0.021 


1623 


n 


12.5 ± 0.6 


1623 


3C249 


n orison 

0.30920 


H 1 


2 


n 


i o f7on i n run 

12.720 ± 0.040 


02 


n 


8.0 ±1.0 


02 


3C249 


0.30809 


H I 


3 


n 


i i ocn i n nor: 

13.860 ± 0.085 


02 


n 


35.7 ± 5.6 


02 


3C249 


0.30788 


H 1 


4 


n 


"i o nr"n 1 n i -in 

13.050 ± 0.140 


02 


n 


,111 1 \ 1 A 1 \ 

22.0 ± 4.0 


02 


3C249 


n on /i <i£? 


H 1 


5 


n 


12.860 ± 0.080 


02 


n 


10.0 ±1.0 


02 


3C249 


0.27021 


H 1 


6 


n 


13.210 ± 0.125 


02 




43 


02 


3C249 


0.26664 


H I 


7 


n 


"i i o on i n n7c 

13.320 ± 0.075 


02 


n 


C7 n I in n 

67.0 ± 19.0 


02 


3C249 


0.26473 


H I 


8 


n 


"i i nnn i n -i nn 

13.900 ± 0.100 


02 


n 


27.5 ± 4.5 


02 


3C249 


0.27021 


Si III 


9 


n 


12.030 ± 0.190 


02 




19 


02 


3C249 


0.24676 


HI 


10 


n 


14.460 ± 0.067 


0203 


n 


37.2 ± 1.2 


0203 


3C249 


0.24396 


HI 


11 


n 


14.200 ± 0.035 


02 


n 


56.9 ±4.7 


02 


3C249 


0.24332 


HI 


12 


n 


13.940 ± 0.080 


02 


n 


52.8 ± 8.7 


02 


3C249 


0.24239 


HI 


13 


n 


13.800 ± 0.055 


02 


n 


30.7 ±3.8 


02 


3C249 


0.23828 


HI 


14 


n 


13.840 ± 0.120 


02 


n 


14.0 ± 1.0 


02 


3C249 


0.24681 


Si III 


15 


n 


12.099 ± 0.125 


0323 


n 


8.4 ±3.2 


0323 


3C249 


0.20093 


HI 


16 


n 


13.250 ± 0.100 


02 


n 


12.0 ±2.0 


02 


3C249 


0.20068 


HI 


17 


n 


13.970 ± 0.090 


02 


n 


38.5 ± 8.9 


02 


3C249 


0.17371 


NV 


18 


n 


13.190 ± 0.130 


02 


n 


24.0 ±7.0 


02 


3C249 


0.19695 


HI 


19 


n 


13.080 ± 0.090 


02 




43 


02 


3C249 


0.18599 


HI 


20 


n 


12.970 ± 0.070 


02 


n 


23.0 ±4.0 


02 


3C249 


0.17593 


HI 


21 


n 


13.740 ± 0.060 


02 


n 


39.0 ±6.0 


02 



NOTE. — Table 5 is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is sho wn he re for 
guidance regarding its form and content. A complete description of the columns in this tabic can be found in Section \2 . 2| of the 



TABLE 6 

Summary of IGM Detections and Results 



Ion 
(1) 


A/total 

(2) 


N 
(3) 


(4) 


(5) 


dAf/dz 3 
(> 10 mA) 
(6) 


dAf/dz 3 
(> 21 mA) 
(7) 


dAf/dz 3 
(> 30 mA) 
(8) 


/9 
(9) 


OVI 


118 


111 


< 0.40 


6.095 


42.0±l° 3 7 


27.7±2;! 


22.2+^4 


2.075 ± 0.119 


NV 


29 


25 


< 0.396 


5.408 


70.41*1;! 


20.6+g-j 




2.001 ± 0.235 


CIV 


36 


29 


< 0.116 


2.523 




14.3111 


u.rtl-1 


1.774 ±0.157 


cm 


51 


50 


< 0.40 


5.644 


23.8+iV 


17.3+^ 




1.859 ±0.092 


Si IV 


34 


30 


< 0.24 


4.317 




11.81" 




1.729 ±0.125 


Si III 


62 


57 


< 0.40 


5.244 




10.4+r* 


74+2.1 


1.739 ±0.102 


HI 


797 


746 


< 0.40 


5.382 


283.9+^:1 


175.2l£° 


144.0±g;g 


1.680 ± 0.030 



Poorly Constrained Ion Species 



FeUI 


15 


11 


< 0.40 


5.853 


CII 


18 


18 


< 0.296 


4.900 


Fell 


5 


1 


< 0.40 


5.913 


Si II 


15 


15 


< 0.372 


5.345 


SII 


3 


3 


< 0.373 


5.352 


SIII b 


1 


1 


< 0.40 


5.924 


SIV b 








< 0.40 


5.889 


NIP 


8 


8 






NIII C 


3 


3 






3 Number 


of absorber 


3 per u 


nit redshift 


path lengt 


b No stat 


sties are pr< 




for this ion 




C Though 


several dctc 


ctions 


of this ion 




a result, s 


11 detection 








therefore r 


efrain from 




ng any such 


statistics 



9 3+ 
+5.3 
2.2 
+ 13.8 
-5.0 
7 1+3.6 
1 - L -1.T 
13.4 
4.8 



9.6 
11.2 



10.0 



5.1 

8.0 
4.3 
7.1 
6.2 



+4.2 
1.8 
+3.7 
1.8 
+ 7.0 
2.5 
+3.6 
1.7 
+ 10.1 
3.6 



3.7 
7.5 
1.3 
6.1 
3.1 



+3.5 
1.4 
+3.5 
1.7 
+7.0 
2.5 
+3.3 
1.5 
+7.1 
2.5 



2.216 ±0.428 
1.587 ±0.121 
1.922 ± 1.505 
1.568 ±0.157 
1.568 ±0.157 



but note th> 



;ed down to the ir 
ick of detections, 
atalog, no system 
derived from tho 
detections for 



ated equivalent ' 



:e dctecti' 



nay be heavily bi; 
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is obtained by correcting f2j on for metallicity and ion- 
ization frac tion. Generally, we assume a metallicity of 
Z = O.IZq (|DS08f ) and an ionization fraction equal to the 
peak collisional ionization equilibrium (CIE) ion abun- 
dance, /ion j listed in Table [31 but we discuss alternate 
choices for the product Z /; on in Section 13.31 The quan- 
tity ^|qm is thus 



(ion) 
IGM 



CPcr 



fim H 



Z(M/H) / io 



1.83 x 10 



-2:i 



dz 
v,2 



-(N ion )d{logN ion )dz 



"70 



Z(M/H) Q / it 




dJV(log N- K 



dz 



(Mon)i 



AlogiV ic 



A;. 



(6) 



where the factor of Azj/z max weights the contribution 

of each bin in dz. Uncertainties in oIq^ are calcu- 
lated as for ri; on ; systematic uncertainties arising from 
our assumptions are not included. Additional systematic 
uncertainties in this quantity likely exist due to effects 
of nonequilibrium ionization, photoionizatio n, or devia- 
tions f rom the assumed 10% solar metallicity. Sh ull et all 
(j 2 1 lh looked into these effects in detail, especially with 
respect to O VI, and we discuss them further in Sec- 
tion urn 

2.4.2. The Lya Forest 

The Lya forest is an important resevoir of baryons, 
with previous surveys finding up to ~ 30% the total 
baryonic content of the universe in this photoiqnized 
phase of the IGM at low redshift (e.g., IPenton et alj2000t 
iLehner et ail 120071; iDanforth fc Shullll2008D . As the cat- 
alog presented here contains the largest sample of Lya 
forest absorbers to date, with a total Lya redshift path- 

( HI) 

length Az — 5.38, it provides useful constraints on S^jqm, 
the fraction of baryons found in the forest. To this end, 
we use thre e methods of c alcula t ing th a t quant i ty: th e 
met hods of IPenton et all (pOOOT l200l . ISchavd {2001), 
and IShull et al.l (1201 lh. 

The method of IPenton etaD (|2000f ) assumes that Lya 
forest absorbers are singular isothermal spheres and fur- 
ther assumes an impact parameter, p, scaled here as 
p = (100 kpc)pioo- As a result, the assumptions of this 
method are only valid for H I absorbers at fairly low 
column density (logiV < 14.5). The method adopts an 
approximate, frequency-integrated H I photoionization 
rate given by, 



Phi ~ (2.49 x 10 



-14 



-V- 



2.-S 



4.6 



(7) 



where J-23 = 1 is the ionizing radiation field at 1 Ry- 
dberg in units of 10~ 23 ergs cm -2 s _1 Hz sr _1 and 
a s is the spectral index of the radiation field, set to 
a s = 1.8. Thus, scaling the H I column density as 



Nm = (10 cm )Nn, we can define 



IGM — 



(1.59 x 10 9 Mq)H (pl 00 J_ 23 [4.8/(a s + 3)]) 



1/2 



cirp 2 p c 



dM(N m ) at1/2mat 
d{N m )dz 



:8.73 x 10" 

E E 

j=0 i=\og N m 



'70 



dz 

J -25 PlOO — 

dN(\ogN m ) 

dz 



4.8 



1/2 



N i/» 



(8) 



'■J 



Our preferred method for low column density ab- 
sorber s, however, is the revised method of IShull et"aLl 
(120111). Althou g h this method retains the assumptions 
of lPenton et~aT1 (2000) regarding the scale and geometry 
of the IGM clouds, it also includes redshift-dependent 
corrections for the evolution in the space density of ab- 
sorbers, 4>{z) oc ( 1 + z) 3 , as well a s the hydrogen pho- 
toionization rate. lHaardt fe Madaul (|2012l ) found a rapid 
increase in the metagalactic ionizing background radia- 
tion from z = to z = 0.7, fitted bv IShull et all (|2011[ ) 
to the form, 

r H i= (2.28 x 10- 14 s- 1 )(l + z) 4 - 4 . (9) 

Thus, scaling the electron temperature as T — 
(10 4 - 3 K)T 4 . 3 , they obtain, 



f^ = (9.oxnr 5 )^ 



-5\ ' l 70 1 Pl00^4 l .3 363 (l + z )°' 2 



[O m (l + z) 3 + fi A ] 1/2 

Nl( 2 d(logN m )dz 



dAf (log A/hi) i /a 



dz 



(9.0 x W-^pltfXf 3 



log N mt 



(l + z,) - 2 



E ^ i , 

j=0 i=logJV min V IPm(l + z j) 3 



CI, 



,1/2 



(W(logN m ) 



dz 



JV 14 _ 



(10) 



Fina lly, we apply the method described by ISchavd 
(|2001D . This formalism assumes that the IGM clouds 
are gravitationally bound and that the observed column 
densities are characteristic over the local Jeans length. 
We expect that these assumptions better characterize 
high column density absorbers (logiVni > 14.5) than do 
the other methods. Hence, scaling the temperature as 
T = (10 4 K)T 4 , they define, 



1/3 rpO.59 



nn = (L46xl0- 4 ) h ^T%Tl 
»" dAf(Nni) 



dz 



Nl( 3 d{N m )dz 



n0.59 

L 4 



(1.46 x 10" 4 ) h^Y^Tl 

Zmax logAf m ax /r 7 » /Vi t 

\p I \dJV {log N m) 



j=0 i=log N„ 



N 



1/3 



14 



(11) 
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where we set T H i = 0.03 x 10 ~ 12 s" 1 at z ~ (|Shull et al.l 
[1991 IWevmann et al.ll200l . 

For all methods, we assume a temperature of T = 
2 x 10 4 K and sum over bins in column density and red- 
shift. We set pioo = 1 and normalize to N14 — 1 through- 
out. For our Az pathlength in Lya, cosmic variance is 
expected to contribute insignificantly to the errors, which 
instead are dominated by Poisson error ()Penton et al.l 

120041) . We therefore ca lculate t he error in ^ the 
one-sided Poisson error (|Gehrelsl[T 986) in dJ\f(d log N, dz) 
bins summed in quadrature. 

2.4.3. Note on Inconsistencies in DS08 

The calculations of contributions to fl by both met- 
als and H I were described and calculated inconsistently 
m DS08 owing to incorrect cosmological factors. The 
baryon fractions listed in their Tables 12 and 13 were 
computed with incorrect cosmological correction s that 
differed from those given in Equations 7-10 of DS08 
and consequently overestimated values of fib for all ions 
discussed. Our new HST archive survey properly in- 
cludes these effects, as well as the redshift evolution of 
the hydrogen photoionization rate, T h(z). The origi- 
nal baryon fraction values tabulated in lDS08l are unsuit- 
able for comparison to other observations or simulations, 
though the absorber statistics that do not depend on dX 
(e.g., dAf/dz) remain valid. Detailed corrections of the 
data tables will be submitted as an erratum. For the sake 
of comparison to the present work, we briefly describe the 
is sue he re. 

IDS08I calculated contributions to fib using a redshift- 
invariant path length dX instead of dz, but the defi- 
nition of dX was applied inconsistently. The text de- 
scribed the conventional "absorp tion pathlength func- 
tion" of lBahcall k Peeblesl (fT969l) : 

dX = (1 + zf [fi m (l + zf + 0\]~ 1/2 dz , (12) 

but the calculations incorrrectly implemented the path- 
length function as 

dX = (1 + z)- 1 [O m (l + zf + n A ]- 1/2 dz . (13) 

The difference in the two formalisms is a factor of (l+z) 3 , 
negligible at z ~ but a factor of ~ 3 at the highest 
redshifts surveyed. As a result, the corrected dN jdX 
values using the proper formalism as presented in the 
erratu m are lower than those published in Table 11 of 
DS08 by 25 - 39% for most ion species. Because C IV is 
only seen at z < 0.1 in the STIS/E140M data, its value 
of dj\f jdX only decreases by roughly 13%. 

The cosmological parameters f2 ion and ^[qm <|DS08l . 
Equations 6 and 8 and Table 12) scale with dJ\f/dX. 
Recalculating the values of Table 12 with the correct dX 
treatment results in a ~ 40% decrease in the values of 
^igm m mos t cases (~ 12% for C IV). Interestingly, this 

~ 40% decrease in ^|qm' ) / fib brings the DS08 value more 
in line with other published values at the same sensitivity 
limits. 

The overall baryon density acco unted for by Lya ab- 
sorbers (Equations 9 and 10 of DS08) suffer similar 

adjustments, changi ng to flj^/f l b = 2 0.7 ± 2.4% if 
calculated with the iPenton et al.l (|2004l ) method and 



fl^J flh = 17.5 ± 2.6% if calculated with the ISchavd 
( 2001) method over a column density range, 12.5 < 
logA^Hi < 16.5. These adjusted values still do not ac- 
count for the redshift evolution of Th(z) oc (1 + z) 4A 
which is included in the method of IShull etahl f|201 ID 
as applied to the present paper's catalog and enters the 
equation as the square root. 

All further re ferenc es in the present paper to values of 
fib reported bv lDS08l refer to the corrected values. 

3. RESULTS AND DISCUSSION 
3.1. Properties of the Catalog 
3.1.1. Comparison to Previous O VI Catalogs 

As previously discussed, a number of studies have been 
undertaken to study the WHIM as traced by FUV O VI 
absorption, most not ably IDS08I . iTripp et all 1)20081 ) , and 
I Thorn fc~C hen (2008j). Each of thse studies broadly used 
the same STIS data sets, albeit with independent reduc- 
tion processes, some minor variation in sight-line selec- 
tion, and with various degrees of assistance from FUSE 
data. Their line selection and measurement method- 
ologies, however, vary considerably, and produce sig- 
nificantly different catalogs. Figure Q] presents column- 
density histograms of the O VI detections reported by 
each of those catalogs, as well as the catalog presented 
here. As these variations can lead to different infer- 
ences and derived quantities, they must be kept in mind 
when interpreting this work or comparing studies of the 
WHIM. Though this discussion focuses on O VI, it is 
analogous to the distributions of other ion species. 

The differences in the catalogs can be attributed pri- 
marily to the line selection criteria as well as approach 
to splitting detected featur es int o subcomponents. In 
terms of detection criteria, DSOl took an inclusive ap- 
proach indexed to Lya, in which any feature detected at 
4 a or gre ater at close to the appropriate wavelength was 
included. DS08's approach to component splitting, how- 
ever, was conservative: most blended lines were reported 
as a single absorption feature and most clearly split metal 
lines were reported together if the associated H I absorp- 
tion was of ambiguous structure. Together, these ap- 
proaches lead to more detections at low column densities 
(some of which may be erroneous identifications), fewer 
detections at intermediate column densities, and a flatter 
histogram when compared to other catalogs. In addition, 
a few O VI absorbers with little or no H I absorption at 
the same redsh i ft wer e missed in that study. In contrast, 
I Thorn fc Chenl (|2008f ) took a very conservative approach 
to inclusion, requiring that both lines of the O VI dou- 
blet be detected blindly without reference to H I. This 
leads to fewer detections overall, especially at low col- 
umn de nsities . Subcom ponents are split m ore liberally 
than in IDS08I . Finally, (Tripp et al.l (|2008[ ) took an ap- 
proach that combined the two aforementioned methods. 
The O VI doublets were identified blindly before a sec- 
ond search pass looked for O VI at redshifts suggested by 
strong absorption from other ions. The issue of subcom- 
ponents was avoided by presenting two sets of measure- 
ments: one that considered absorption "systems" with 
no splitting and one that split the systems into all dis- 
cernable components. As a result, the histogram for their 
"component" detections is strongly peaked at intermedi- 
ate column densities, while their systemic measurements 
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Fig. 1. — Histograms comparing O VI detections as reported in this work, [DS08, Trip p et all ^2008fl. and IThom k, C hen (2008). Bins 
in logTVovi are 0.2 dex in width, centered on the plotted points. The two separate sets of points for lTrip p et al.l 112008!) re present their 
"system" (Table 2) and "component" (Table 3) measurements. Our catalog and component results of ITripp et ail 1 12008!) are strongly 
peaked at logJVovi ~ 13.6. 



lead to a distribution shifted to high column densities 
and reduced in total number of detections. 

As described in Section 12. 2\ our catalog is interme- 
diate in its approach to component splitting compared 
to the aformentioned studies. It includes all reason- 
able detections included in previous catalogs that are 
not contradicted by new COS data. As a result, it is 
the largest catalog in terms of number of H I and O VI 
detections. The O VI distribution is strongly peaked at 
logiV ks 13.2, and it contains an intermediate num ber of 
det ections at lower co lumn densities compa red to IDS08I 
and lTripp etaLl (|2008H . Compared to lDS08l it corrects a 
number of misidentifications at low column densities and 
augments the number of detections in light of new data 
and work by other research groups. Our catalog has 96 
more H I absorbers and 28 more O VI absorbers than 
DS08. The increase can be attributed to the inclusion of 
new absorbers and sight lines as well as the more liberal 
splitting of absorbers into separate features. The cor- 
rected misidentifications were primarily weak metal lines 
mistaken for other species (e.g., Si III mistaken for Lya) 
or uncommonly seen metals associated with high column 
density H I absorbers mistaken for H I. As a result, this 
catalog supercedes that of [PS08 for most purposes. 

3.1.2. Contamination in reported absorbers in STIS data 

Our catalog contains 381 H I detections in the 20 sight- 
lines for which both COS and STIS and/or FUSE data 
are present. Of these detections, only 330 are detected 
by COS as well as STIS, as depicted in Figure[2j In most 



STIS+COS Targets 
STIS + COS Targets, no contaminants 





log N(HI 



Fig. 2. — Histogram comparing detected H I absorbers for sight 
lines with both STIS and COS data. The solid line represents all 
detections, while the dotted line represents detections seen in both 
data sets. 

cases, because the COS data have much higher S/N than 
the STIS data, this discrepancy suggests that some ab- 
sorbers previously reported in the literature are in fact 
noise features rather than IGM absorption lines. For col- 
umn densities N m < 10 132 cm~ 2 , between 20% and 40% 
of STIS detections are not supported by the COS data. 
In the column density range 10 132 — 10 136 cm -2 , this im- 
proves to around 5% to 10% contamination, while higher 
column densities are nearly entirely free from contami- 
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nation. These results suggest that censuses of the low- 
rcdshift Lya forest at column densities below 10 13 5 cm" 2 
may suffer significant contamination in the form of line 
misidentification that is not accounted for in standard 
completeness corrections. Such contamination under- 
scores the need for high-S/N observations of the IGM 
sensitive to Lya equivalent widths W\ ~ 5 mA, corre- 
sponding to logA^Hi ~ 12.0. 

3.2. Properties of the Absorbers 

Table [5] summarizes our detection statistics for several 
ion species, listed in Column (1). The ions in the bot- 
tom half of the table have poor statistics, and should be 
viewed with skepticism. Columns (2) and (3) report the 
total number of 4 a detections in STIS data and the num- 
ber of detections with false positives indicated by COS 
data removed, respectively. The latter number is the set 
of detections used in all statistics and derived quantities. 
Column (4) lists the range of redshifts that the survey 
covers for a given ion, while Column (5) lists the maxi- 
mum redshift pathlength sampled at any sensitivity for 
the ion. Columns (6) through (8) report the absorber 
frequency per unit redshift pathlength integrated down 
to three equivalent widths. The survey is greater than 
50% complete at the 30 mA scale, but at lower equiva- 
lent widths the completeness corrections take on a larger 
role. Finally, Column (9) provides a fit, in terms of the 
index 0, to the distribution of absorbers as a function of 
column density per unit redshift, 



d 2 N(> N) 
dz dlogN 



oc N~ 



(14) 



where we have used an error-weighted fit to the cumula- 
tive distributed histogram of dN/dz. This distributed 
approach may systematically flatten the sl ope, , by 
0.05 — 0.1 dex relative to simple histograms (DSOl). 
In Table [H we summarize our results for the baryon 

fraction for metal ions, ^iqm an d ^ion, listed for equiva- 
lent widths integrated down to 10 mA and 30 mA. 
is scaled by metallicity and peak ionization fraction as 
described above, and it is reported relat ive to the total 
baryo n fraction, fl b = 0.0455 ± 0.0028 (jKomatsu et all 

lm|). 

3.2.1. Individual Species 

HI. — Our results for the column-density distribution 
of H I absorbers are presented in Figure [3] The his- 
togram and dN/dz plots both exhi bit a d ip around 
logA^Hi ~ 13.5, which was also seen in IDS081 The most 
likely explanation for this feature is that it is an arti- 
fact of the measurements arising from the detectability 
of higher-order Lyman lines. There is therefore more fre- 
quent use of COG measurements at column densi t ies be - 
yond this point. As discussed in iDanforth et al.l (2006), 
measurements of absorbers that exhibit only Lya may be 
biased toward lower column densities. This effect may 
cause some absorbers to scatter into lower bins, creating 
the dip in the distribution near the transition between 
Lya-only measurements and COG measurements. We 
expect that this effect has little impact on our absorber 
statistics, since it does not affect the total number of 
detections. A second dip is seen in the dN '/dz plot at 



logTVjji ~ 12.8. It is not significant at the 95% confi- 
dence level, and the incompleteness correction dominates 
in this regime. 

The entire distribution yields a power-law ind e x, = 
1.68 ± 0.03, in agreement with IPenton et all (120041) , 
who found = 1.65 ± 0.07; ILehner et all ([20071 ). who 
found = 1.76 ± 0.06; and DS081 . who found = 
1.73 ± 0.04. Limiting the fit to bins that are greater 
than 50% complete, however, yields a slightly steeper 
value, = 1.74 ± 0.04, which may indicate system- 
atic effects from our completeness corrections. The 
completeness-corrected line frequency integrated down 
to 30 mA (log Nm sa 12.74) is dN/ d z = U 4± 7 6 , in good 
agreement with lPenton et all (|2004f ). IDS08I found a lower 
value of dN/dz = 1291^, where the discrepancy likely 
arises from the methodological differences described in 
Section EXlJ 

The distribution of the Doppler parameter, b, for 
H I provides constraints on the thermal state of 
the IGM. The distribution of measured 6-values, lim- 
ited to absorbers with errors in b and Nm less than 
50%, is shown in Figure 21 with a range from b = 
6 — 137 km s _1 with a mean, median, and standard de- 
viation of 34, 30, and 17kms _1 , respectively. Values 
of b = (30 km s _1 )&3o correspond to a thermal tempera- 
ture of T w (54, 500 K)6| , but this number can be taken 
only as an upper limit because turbulence, bulk flows, 
and undetected subcomponents will inflate the measure- 
ments. The broad Lya absorbers (BLAs) predicted by 
simulations and oft en associated with the WHIM (e.g., 
IRichter et al.l I2006D are difficult to characterize in the 
noisy STIS echelle data. Roughly 25% of our H I cat- 
alog has a measured b > 40 km s _1 , but the errors are 
large in this regime. In addition to the large statistical 
errors, systematic effects such as unresolved subcompo- 
nents, insufficiently characterized AGN continua, and in- 
strumental artifacts likely add further uncertainty to the 
upper end of the distribution of b- values. The treatment 
of BL As i s dealt with more car efully in ILehner et al.l 
(|2007f) and IDanforth et all (|2010f ). with the latter utiliz- 
ing much higher S/N data from COS to further constrain 
their properties; we refer the interested readers to those 
works. 



O VI. — We find a power-law index for the entire 
O VI distribution (Figure top panels) of = 2.08 ± 
0.12, in good agreement with o ur group's previous work 
(|Danforth & Sh ull 2008, 120051) . Interestingly, we see 
a turnover in the dN/dz column density distribution 
at logVovT ~ 1 3.5. This effect was also seen by 
Danforth & Shull (2005;), albeit with just five abs orbers 
at such low column densities. On the other hand, IDS08I 
found no evidence of a turnover, with 26 absorbers in the 
bins at log Vovi < 13-5. We report 29 absorbers in our 
three lowest bins (12.9 < logiVovi < 13.5) over a larger 
pathlength in z. The turnover seems to arise from the 
combined effect of a stronger peak in the distribution of 
absorbers at log Vovi ~ 13.6 coupled with an increase 
of onl y 3 absorbers with logVovi < 13-5 compa red to 
IDS08I . Our correction of several absorbers from the lDS08l 
catalog that were misidentified as O VI explains the lack 
of a large increase in the number of absorbers in these 
bins despite our larger pathlength. Despite its statistical 
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Fig. 3. — Distributed histograms of H I detection statistics and A/"(HI) = dj\[ /dz. Left panel shows the number of absorbers per 0.2 dex 
bin, uncorrected for completeness. Right panel shows the completeness-corrected differential dj\f /dz vs. log JVhi. Error bars in both panels 
include one-sided Poisson errors, while the right panel errors also include contributions from Az uncertainty. 
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Fig. 4. — Distribution of Doppler 6-parameters of 626 of 746 H I absorbers, limited to absorbers with errors in b and ./Vhi less than 50%. 
The median and sample standard deviation of the distribution are 30 km s" 1 and 17 km s — , respectively. Because some absorbers may 
actually be multiple blended components, this distribution provides upper limits on b-values. 



significance, the turnover should be treated with cau- 
tion. The statistical errors on individual measurements 
at these column densities are quite large, and the method 
of catalog construction may introduce systematic biases 
at low column densitie s owin g to the predominant re- 
liance on the catalog of IDS08I in that regime. It is also 
worth noting that our completeness corrections exceed 
50% in these bins. A new COS survey at higher S/N 
will be required to conclusively address the issue of this 
turnover. If the turnover is real, it may be more realis- 
tic to apply the power-law formalism only to data with 
logA^ovi > 13.5; this approach yields a much steeper 
slope, $ = 2.50 ±0.26. 

Our results for the frequenc y of O VI absorber s com- 
pare favorably to simulations. [Smith et al.l (|201 If ) inves- 



tigated the distribution using the N-body hydrodynamic 
code Enzo in a 50/i -1 Mpc comoving box with 1024 3 grid 
cells (run 50_1024_2 in that work). Cooling in CIE and 
with both collisional and photoionization (C+P) were in- 
vestigated. We plot their results for the cumulative dis- 
tribution of dJ\f(> N(yvi)/dz against our empirical dis- 
tribution in Figure |5] Though CIE adequately explains 
the distribution at high column densities, it over-predicts 
the distribution at lower column densities, suggesting an 
increasing role of photoionization at these scales. These 
smaller clouds may reside in areas of the IGM that are 
characterized by overall lower baryonic densities, mak- 
ing them more sensitive to the effects of photoionization. 
Alternately, perhaps they reside in the circumgalactic 
medium, where higher UV radiation helds may increase 
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Fig. 5. — Simple histograms of detection statistics and distributed histograms of dAf /dz vs. logAf for O VI (top), C IV (middle), and 
N V (bottom). 



the role of photoionization. 

In Figure [3 we present the O VI Doppler 6-parameter 
distribution, again limited to absorbers with errors in 
b and iVovi less than 50%. Nearly all measurements 
in this distribution were performed using the the Voigt 
profile routine described previously, which accounts for 
instrumental resolution. The distribution spans a range 
of b = 6 — 78 km s _1 with a mean, median, and stan- 
dard deviation of 32, 28, and 16 km s _1 , respectively, in 
excellent agreement with IDS08I . Q ualitatively, t he dis - 
tribution is similar to Figure 13 of lTripp et all (l2008h . 



although they find a higher density of absorbers with 
b < 10 km s -1 . We suspect that this discrepancy arises 
from their more aggressive splitting of systems into sub- 
components, which may artificially lower their distribu- 
tion of 6-values. Eleven absorbers (14%) in our catalog 
have 6-values corresponding to temperatures at or below 
10 5 - 45 K, the peak in CIE. Only two (3%) correspond to 
temperatures less than 10 5 K as expected from photoion- 
ization. 

To investigate the correlation of O VI absorption with 
H I absorption, we combined all components of each ion 
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Fig. 6. — Cumulative dj\f{> N)/dz for O VI absorbers com pared 
to simulation results of run 50_1024_2 from lSmith et al.1 |(20Tj) fea- 
turing collisional and photoionization (C+P - solid line) and CIE 
cooling (dashed line). 

species with |Av| < 50kms _1 , yielding "systemic" col- 
umn densities and velocity displacements. We then cor- 
related the two set ol absorbers, counting any O VI ab- 
sorption with a redshift within several Av offsets of the 
rcdshift of the H I as being associated with that H I 
absorber. The resulting distributions are plotted as his- 
tograms as a function of column density in Figure [8] The 
strongest H I absorption systems nearly always have as- 
sociated O VI. In all but a handful of cases, O VI has 
associated H I absorption. We found four cases of O VI 
absorption without apparent associated H I absorption 
within \Av\ < 100 km s _1 and only one case without 
H I absorption within \Av\ < 300 km s -1 . This kine- 
matic association confirms that H I and O VI often arise 
in spatially asso ciated, multiphase sys t ems, as predicted 
by simulations (iCen fc Ostrikerl [19991: iDave et all 1 2001 1 : 
Smith et al]l20 1 11) and observed jDanforth fe Sliuillbooa 
Tripp et al.ll200l iSavage et alll2011bl) . 

C IV and N V. — The Li-like C IV and N V ions have 
also been suggested as possible tracers of WHIM gas, 
although photoionization probably also plays a role for 
these species. We report 29 C IV absorbers, shown 
in Figure [5] (middle row), but our survey is sensitive 
to C IV over a much shorter total redshift pathlength 
(Az = 2.52) compared to other ions. The distribution 
yields a power-law index of (3 — 1.77 ± 0.16 and an ab- 
sorber frequency integrated down to 30 mA (A1548.2; 
log N C w ~ 12.87) of d M/dz = 12±£ Our C IV catalog 
has changed little since IDS08I . so, unsurprisingly, these 
results are in good agreement with the results from that 
study and its discussion of the ion. As in the O VI distri- 
bution, we see a turnover in the dM /dz column density 
distribution of C IV at logiVciv ~ 13.2, but it is not 
significant at the 95% confidence level. We invoke all of 
the same caveats in ascribing meaning to this feature as 



we did for the O VI results. 

Despite having more than double the redshift path- 
length of C IV, the N V ion has only 25 absorbers, 
shown in Figure [5] (bottom row) , yielding a power- law 
index of /3 = 2.00 ± 0.24 and an absorber frequency in- 
tegrated down to 30 mA (A1238.2; logiV NV « 13.15) of 
dN /dz — IOI3. The low frequency of detections likely 
is attributable to the low abundance of nitrogen com- 
pared to oxygen or carbon. As a result, most detections 
are weak, and it will remain difficult to draw many con- 
clusions from N V absorption until a systematic survey 
in more sensitive COS data is undertaken. The slight 
turnover apparent in the dM /dz column density distri- 
bution is not statistically significant. 

C III, Si III, Si IV, Fe III and other ions. ~ C III and Si III 
are expected to arise in photoionized conditions, and they 
are well characterized by the high-resolution STIS data, 
owing to their relatively narrow line widths and strong 
absorption f eatures . Our catalog's primary differences 
from that of DS08 arise from a more careful treatment 
of the narrow subcomponents seen in many of these ab- 
sorbers. We report 50 C III absorbers, shown in Figured 
(top row), yielding a power-law index of — 1.86 ± 0.09 
and an absorber frequency integrated down to 30 mA 
(logAcin « 12.67) of dM/dz = UA±% 7 2 . We find 57 
Si III absorbers, shown in Figure^ (second row), yield- 
ing a power-law index of (3 — 1.74±0.10 and an absorber 
frequency integrated down to 30 mA (logiVsan « 12.14) 
of dM/dz = 7A±f l, We see no significant turnover in 
their dM /dz column density distributions. 

We report 30 Si IV detections, shown in Figure ® (third 
row), which follow a distribution quite similar to that 
found in I)S(Ks. They found a power-law index of /3 = 
1.92 ± 0.17 that was somewhat steeper than that of the 
other low-ionization state species. Our catalog yields an 
index of (3 — 1.73 ±0.13, bringing it more in line with 
the other ions. We find an absorber frequency integrated 
down to 30 mA (A1393.8; log A si i V » 12.53) of dM /dz = 

9 2+ 45 
1 g 'i 

IDS08I made an ambitious attempt to characterize Fe III 
absorption in the IGM, reporting 14 detections compris- 
ing a steep distibution with j3 — 2.2 ±0.4. However, they 
expressed some caution about the validity of these re- 
sults, owing to the high Fe III abundance ratios. Our re- 
analysis of the data suggests that all Fe III results should 
be treated with skepticism. Although we report 11 de- 
tections at the nominal 4 a level, nearly all detections for 
which COS data were available were revealed as misinter- 
preted noise features. It is thus likely that a significant 
number of those 11 absorbers may vanish, too, if ob- 
served at a higher S/N with COS. In reporting our Fe III 
statistics for this species in Table HI we tentatively note 
that they are poorly constrained. Table [5] further reports 
our detections of several other low-ionization species. As 
with Fe III, we mark them as poorly constrained owing to 
the paucity of detections. No statistics are presented for 
N II (A1084.0, A915.6) or N III (A989.8) because no sys- 
tematic search was made for these species. The listed de- 
tections were primarily serendipitous corrections of past 
misidentifications in the literature, and they are primar- 
ily associated with high column density H I absorption 
systems. 
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Fig. 7. — Distribution of Doppler 6-parameters of 79 of 111 O VI absorbers, limited to absorbers with errors in b and N less than 50%. 
The median and sample standard deviation of the distribution arc 28 km s —1 and 16 km s , respectively. Because some absorbers may 
actually be multiple blended components, this distribution provides upper limits on b-valucs. 

variations in the data itself. That work also suggested 
that perhaps their values (9%, 5%, and 3% for O VI, 
N V, and C IV, respectively) represented a true gradient 
the WHIM baryon density as a function of temperature, 
because their derived densities follow the trend in peak 
CIE temperatures. This interpretation is not supported 
by our improved values, though we are unable to rule it 
out. 

The largest systematic uncertainties in these values 
arise from the metallicity and ionization fraction assump- 
tions. In the above values, we assume a metallicity of 
10% solar and an ionization fraction corresponding to 
peak CIE value. Recent simulations (jShull et al.ll2011[) 
suggest that the weighted mean product of these quan- 
tities may be closer to fovi(Z/Z@) — 0.01, half of our 
assumed value of 0.02. Using this number in our cal- 
culations effectively doubles the baryonic content traced 

by O VI, yielding ~ 17 - 2 ± 1A %- IShull et al.l 

(|2011l) also report a column density dependent fit to this 
quantity: fovi(Z/Z Q ) = (0.015) [iV OV i/10 w cm- 2 ] - 70 . 
Evaluating that fit separately for each bin in our 
calculations raises the baryonic content still further, to 
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Fig. 8. — Histograms of all "systemic" H I absorbers (see Sec- 
tion |3]2TL0] for details), in which subcomponents within 50 km s _1 
of each other have been combined (solid line), and systems 
with O VI absorbers within 50 km s — 1 (dotted line), 100 km s — 1 
(dashed line), and 200 km s — 1 (dash dotted line) of the systemic 
H I redshift. 

3.3. The Baryon Density of the IGM 

Using the assumptions for metallicity and ionization 
fraction described above, the three best tracers of the 
WHIM yield estimates of the fractional baryon content 

that they trace as re 9 ± 1% (O VI), 6 ± 3% 

(N V), and 7 ± 1% (C IvV lDS08l offered several interpre- 
tations for variations in these values. Most obviously, the 
differences may arise due to dominant systematic errors 
arising from our assumed CIE ion fractions or metal- 
licities, or they may arise simply due to line sensitivity 



IGM 



/n b w 18.9 ± 2.0%, but the scatter in the foviZ 
fit is roughly 0.35 dex. Simulations suggest, however, 
that the actual contribution of O VI to f2 , as calculated 
by the direct summation of O VI within the simulation, 
may be as much as 10 — 40% higher than that derived 
from absorption-line studies, as traced by the creation 
of sy nthetic spectra from the simulation ([Smith et al.l 
120111 ) . We can apply the analogous results fr om the same 
simula tion, derived from run 50_1024_2 from lSmith et al.1 
(|2011|) (B.D. Smith, private communication), to our N V 
and C IV results. Using the column density fits for these 
ions, we obtain baryon fractions of 10.1 ± 6.0% (N V) 
and 22.3 ± 3.7% (C IV). The C IV results change most 
dramatically, possibly because it is likely to be more sen- 
sitive to photoionization effects. 

The baryon content of the ph otoionized Lya fores t has 
been well studied elsewhere. iPenton et al.l (|2004| ) ac- 
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Fig. 9. — Simple histograms of detection statistics and distributed histograms of dAf /dz vs. logiV for C III (top), Si III (second row), 
Si IV (third row), and N V (bottom). 



counted for roughly 29% of the loc a l bary on mass down to 
logTVni = 12.5, and lLehner et al.l (|2007t ) found a similar 
estimate of around 30%. We present our own calcul a- 
tions based o n the metho d s desc ribed in Section [2.4.21 in 
Table The iShull et al.l (|2011[ ) method ismost appro- 
priate at lower column densities, while the lSchavel (2001) 
method, which assumes gravitationally bound clouds, 
is likely more appropriate at higher column densities. 
In estimating a total ratio ^ ya /ftb, a cutoff between 
the two methods should be assumed. Taking a cut- 
off of logA^Hi = 14.5, corresponding to the transition 



in the slope of the column density distribution seen by 
IPenton et all (pOOl . we find fi Ly aM = 20.9 ± 2.0%. If 
the cutoff is instead drawn at logiVHi = 16.5, we find 
= 23.7 ± 2.2%. Because this calculation is 
model dependent, systematic errors likely dominate over 
the statistical errors that we report. 

The baryon fractions reported in Table [8] for the lower- 
ionization species likely trace a metal-enriched subset 
of the gas in the Lya forest. As a result, they are 
probably photoionized, rendering assumptions regarding 
the ionization fraction problematic. Hence, the values 
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TABLE 7 

Baryon Content of the Local Lyc* Forest 



log Nm Range 


M 




OLyc/fii, (%) a 




OLya/Oi, (%) b 


Or c 


^LyaM (%) C 


12.0- 


12.5 


20 


0.0055 ± 0.0043 


12.0 ±9.5 


0.0038 ± 0.0041 


8.3 ±9.0 


0.0036 ± 0.0041 


8.0 ±8.9 


12.5 - 


13.5 


497 


0.0044 ± 0.0010 


9.8 ± 2.2 


0.0039 ± 0.0008 


o.D ztz i.o 


u.uuoo zn u.uuuo 


Q A _|_ 1 *7 
0.^± ZXZ -L. 1 


13.5 - 


14.5 


244 


0.0031 ± 0.0003 


6.7 ±0.6 


0.0038 ± 0.0003 


8.4 ±0.8 


0.0038 ± 0.0003 


8.3 ± 0.7 


14.5 - 


15.5 


49 


0.0012 ± 0.0002 


2.7 ±0.5 


0.0023 ± 0.0004 


5.0 ±1.0 


0.0022 ± 0.0004 


4.9 ±0.9 


15.5 - 


16.5 


3 


0.0001 ± 0.0001 


0.3 ±0.2 


0.0003 ± 0.0002 


0.7 ±0.4 


0.0003 ± 0.0002 


0.7 ±0.4 


12.5 - 


14.5 


671 


0.0075 ± 0.0010 


16.5 ± 2.3 


0.0077 ± 0.0009 


17.0 ±1.9 


0.0076 ± 0.0008 


16.6 ± 1.8 


14.5 - 


16.5 


52 


0.0013 ± 0.0002 


2.9 ±0.5 


0.0026 ± 0.0005 


5.7± 1.0 


0.0026 ± 0.0005 


5.7 ± 1.0 


16.5 - 


19.0 


3 


0.0006 ± 0.0004 


1.4 ±0.8 


0.0035 ± 0.0022 


7.8 ±4.8 


0.0035 ± 0.0021 


7.7 ±4.7 


12.5 - 


16.5 


723 


0.0088 ±0.0011 


19.4 ±2.3 


0.0103 ±0.0010 


22.7 ±2.2 


0.0102 ± 0.0010 


22.3 ± 2.1 


12.5 - 


19.0 


726 


Shull method for log JVhi < 16.5, Schaye method for 


log N m > 16.5: 


0.0108 ±0.0010 


23.7 ±2.2 



a Method of lPenton et al.l 120001) 

b Method of l^jchavd ^OQjj ) 

C Method of lijhull et al.l <2ulH ) 



TABLE 8 

Summary of IGM Baryon Fraction Results 
Ion n ion (10- 8 ) ^gmAV 

(> 10 mA) (> 30 mA) (> 10 mA) (> 30 mA) 

OVI 48.22 ±3.77 40.33 ± 4.28 0.086 ± 0.007 0.072 ± 0.008 

NV 4.19 ± 2.25 2.42 ± 2.85 0.059 ± 0.032 0.034 ± 0.040 

CIV 20.11 ±2.01 19.83 ±2.06 0.067 ± 0.007 0.066 ± 0.007 

CIII 10.83 ± 1.38 8.75 ± 1.46 0.013 ± 0.002 0.010 ± 0.002 

Si IV 11.07 ± 1.57 10.78 ± 1.78 0.099 ± 0.014 0.097 ±0.016 

Si III 17.53 ±3.21 17.26 ±4.27 0.061 ±0.011 0.060 ± 0.015 

Felll 40.16 ± 14.31 29.48 ± 16.93 0.082 ± 0.029 0.060 ± 0.035 

a Sca led by / ion , Z, and Hp; tt b = 0.0455±0.0028 (Ko matsu et al.l 
1201 If ). See text for details. 



of f2jQ^j/J7b for these ions are extremely uncertain be- 
yond the statistical errors given in the table. Recalcu- 
lating them with co lumn density dependent fits from the 
ISmith et all (|201lD simulation gives values 2 to 3 times 
higher, but this huge variation likely indicates imperfec- 
tions in the combined ionization calculations in simula- 
tions to date. 

4. CONCLUSIONS AND SUMMARY 

We have presented a comprehensive, critically evalu- 
ated catalog of low-redshift IGM absorbers present in 
archival legacy STIS and FUSE data. New data from 
COS informed the evaluation, revealing past errors in the 
interpretation of absorption lines. We analyzed 44 sight- 
lines, yielding 746 H I absorbers over a total Lya redshift 
pathlength of Az — 5.38. In addition, we reported on 111 
O VI absorbers, 29 C IV absorbers, and made numerous 
detections of other ion species (N V, C III, Si III, Si IV). 

Our O VI results found distributions comparable 
to those found in past work. We found a dM/dz 
tu rnover at low column de nsities similar to th at seen 
by IDanforth fe"Slnilil (pOOl but absent in iPSOl Using 
a more sophisticated baryon accounting technique than 
previous work, we estimated from our empirical distri- 
bution of absorbers that the WHIM contributes roughly 
19% of the total baryonic density, fif,. We discussed sev- 
eral significant systematic uncertainties in deriving such 
a quantity, most notably corrections for the metallicity 



and ionization state of the IGM. We presented similar 
distributions and quanitities derived from N V and C IV 
detections, of which we reported 25 and 29 absorbers, 
respectively. Similar statistics were reported for other 
ion species at lower ionization states, and we summa- 
rized these results in Table El In particular, we found 
that past studies of Fe III absorption in the low-redshift 
IGM suffer from significant contamination due to line 
misidentification. 

Our catalog of 746 H I absorbers follows a similar dis- 
tribution to that of our group's past work in this field. 
IDS08I used an incorrect cosmological factor to derive 
baryon densities from the detections; we corrected that 
error with our new catalog, accounting for roughly 24% 
of Hi, with Lya forest detections. Taken together with 
inferred WHIM baryon densities, our catalog accounts 
for ~ 43% of in the low-redshift IGM. Unfortunately, 
the present catalog is unable to adequately probe broad 
Lya absorbers, which are expected to trace significant 
baryonic content in the WHIM. Such absorbers require 
high S/N data, better supplied by future work with COS. 

This paper has analyzed most of the IGM absorption- 
line systems at z < 0.4 observed by STIS. With COS, one 
can take these studies much further through a spectro- 
scopic survey of the structure, metallicity, and physical 
conditions of the IGM and CGM. Observing the relation 
of IGM absorption systems to galaxy distributions will be 
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achieved through uniformly processed, high-S/N spectra 
and an unbiased sample of IGM and CGM absorbers, in- 
cluding the Lya forest, metal-line systems, Lyman-limit 
systems, and partial-limit systems. We have proposed a 
sensitive (S/N ~ 30) HST Spectroscopic Legacy Survey 
of the low-redshift IGM and CGM toward AGN selected 
for UV brightness and total redshift pathlength. These 
absorbers will be compared with deep galaxy surveys, 
thereby connecting galaxies to large-scale structure and 
inferring the extent of feedback. These projects will pro- 
duce a complete legacy survey of gas and (L > Lsmc) 
galaxies, located within 1-2 Mpc of the (z = 0.1 — 0.4) 
absorbers. The high UV sensitivity of COS provides an 
excellent tool for studying gas in multiple thermal phases 
and metallicities in the low-redshift IGM and CGM. 
These studies are best done at low redshift (z < 0.4), 
allowing us to use correlations of the absorbers with dis- 
tributions of galaxies, particularly those of low mass and 
luminosity (L < L*) which contribute significant radia- 
tive, thermal, and chemical feedback. While the catalog 
presented here represents an improvement to our census 
of IGM baryonic content, these future studies with COS 



will be able to address conclusively the distributions of 
absorbers at low column densities and broad line widths. 
Poorly constrained ions such as Fe III and N V may be 
probed with more reliability, and the increase in path- 
length will dramatically improve the statistics of all ion 
species addressed here. 
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